Inflation as a Solution to the Early Universe Entropy Problem? 
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There exists the 'entropy problem' of the early universe, that is, why did the universe begin with 
an extremely low entropy and how did it evolve into such high entropy at late times? It has been 
long believed that inflation cannot be the solution since it requires an extremely low entropy to ever 
occur. However, we point out that since the inflation is always accompanied with a horizon, the 
correct probability of inflation is associated with the quantum entanglement entropy, which should 
in principle be larger than what considered previously. This motivates us to reexamine the issue by 
computing the evolution of the cosmological entanglement entropy in the early universe. We invoke 
a toy model of nonlinear generalized Chaplygin gas (GCG), which has the advantage of providing 
a smooth and unitary transition between the inflation epoch and the radiation dominant era. We 
found that soon after the onset of the inflation, the total entanglement entropy rapidly decreases 
to a minimum, and it rises monotonically afterwards throughout the remainder of the inflation and 
the radiation epochs. This indicates that the universe does not need to begin with an extremely 
low entropy; its smallness can be naturally induced by the dynamics of inflation itself. We believe 
that our computation largely captures the essential feature of entropy evolution and can provide us 
insights beyond the toy model. 

PACS numbers: 



I. INTRODUCTION 

An increasing amount of observational evidence sup- 
ports the notion that the early universe has undergone 
an epoch of inflation. However we are still far from under- 
standing the underlying assumptions and resolving some 
of its most crucial issues. One of the open questions 
about inflation is the entropy problem [1] [2] [3] [4] : why did 
the universe start from an extremely low entropy initial 
state that in turn can later give rise to the arrow of time 
in accordance with the second law of thermodynamics? 

Thermodynamic arrow tells us that a system will 
evolve towards the increase of entropy, which simply 
states that a system is more likely to evolve from a state 
with low occurrence probability to that with a higher 
probability. Since the universe is old, in order to ex- 
plain the thermodynamic arrow of time the entropy in the 
early universe must therefore be extremely small, which 
requires a special initial condition. 

On the other hand, the standard inflation paradigm as- 
sumes the initial condition of the universe to be realized 
generically. This conflict [4] becomes more acute if we re- 
gard the initial entropy as the probability for spontaneous 
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formation of a homogeneous domain in the inflationary 
universe. An extremely low initial entropy required by 
the thermodynamic arrow then suggests that the initial 
universe only has a very small probability to evolve into 
the current universe through inflation. 

Several solutions have been suggested in the litera- 
ture. One approach is the spontaneous eternal inflation 
[3][S], where inflation is assumed to occur both forward 
and backward in time. Another is the so-called bubble 
cosmology [6|, where universes were induced as a realiza- 
tion of string landscape[7][8][9]. 

In this paper, we seek a solution to the entropy prob- 
lem within the present framework of the standard model 
of cosmology without relying on assumptions like sym- 
metric time arrows or string theory. The framework will 
only be effective, in that new physics is expected to oc- 
cur near the initial singularity predicted by General Rel- 
ativity, and detailed mechanism for the transition from 
inflation to a radiation-dominated universe is still not 
well- understood. 

Our aim is to show that a low entropy state can be 
dynamically produced from inflation itself, followed by 
the generation of a large amount of entropy. This pos- 
sibility has long been refuted in the previous literature 
[TU] [TTj [T5] based on the argument that inflation itself 
requires an extremely special initial condition to occur. 
However, we hold the point of view that whether the ini- 
tial condition is special or not depends largely on how 
one measures the probability of the degrees of freedom. 
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For example, in the recent development of loop quantum 
cosmology 14] it was found rather differently that the 
probability for inflation to happen is much higher using 
their measure. We should like to point out an important 
fact, which was overlooked by those previous authors, 
that inflation is accompanied with a horizon. Therefore 
the correct measure of the probability for inflation should 
be given by quantum entanglement entropy. Since the 
presence of the horizon increases the entanglement en- 
tropy by the entanglement between the observable and 
the unobservable regions, inflation should have a higher 
probability to occur. This insight motivates us to seek 
new dynamical solution to the entropy problem within 
the inflation paradigm. 

We notice that in the literature there are examples of 
decreasing entanglement entropy [E][16] due to the sit- 
uation where the system is out of equilibrium or simply 
from long-range quantum correlation. In particular, a 
decrease of tripartite entropy after inflation is shown in 
[16] , which inspires us to consider the evolution of bipar- 
tite entropy from quantum entanglement. Indeed, in a 
recent effort by two of the present authors (PC & YN) 
|17j . it was found that a minimum of the entanglement 
entropy induced from the entanglement between cosmo- 
logical perturbations is reached at an early stage of the 
inflation. The calculation made there, however, was clas- 
sical. Since we expect quantum effects to be prominent 
in the beginning of the universe, each degree of freedom 
must subject to quantum fluctuations. We therefore re- 
address the issue properly in this paper. 

This paper is organized as follows. In section |TT] wc 
briefly summarize the reasons and methods for consider- 
ing entanglement entropy as entropy budget of early uni- 
verse, which consist of homogeneous and inhomogeneous 
parts. In section |III| we introduce the effective model 
used in our calculation. Detailed treatment and numeri- 
cal results are summarized in next two sections IV andfV] 
which together consist of the total contribution to entan- 
glement entropy. To make comparison with the previous 
result in |17j . the computation therein is reproduced in 
section [V] for completeness. Results in each section are 
discussed, and the main conclusion is summarized in the 
last section. 



II. ENTANGLEMENT ENTROPY AS 
ENTROPY BUDGET 

In the present framework, most conventional consider- 
ations of entropy are thermal, statistical and entangle- 
ment entropy accounted by von Neumann. Since early 
universe experienced accelerated expansion, it is far out 
of equilibrium, and thermal entropy is ill-defined for lack 
of a universal temperature. Statistical mechanics is also 
not directly applicable on present situation due to long- 
range nature of gravitational force. 

On the other hand, von Neumann entropy from 
quantum entanglement is perfectly well-defined in non- 



equilibrium system, and can cope with problems statisti- 
cal entropy encounters in gravitational systems like loss 
of information. Moreover, von Neumann entropy fits 
into unitary framework of entire universe in that it mea- 
sures information encoded in entanglement between sub- 
systems, while thermal entropy is of dissipation nature. 
Since the early universe can be described as a quantum 
system without strong decoherence, it is justifiable to in- 
voke von Neumann entropy to find the evolution of total 
entropy of the universe. We will thus use the quantum 
entanglement entropy instead of thermal entropy in the 
following discussion to help clarify our setup and to quan- 
tify the amount of information in a quantum system. 

Here we will clarify what we mean by von Neumann en- 
tanglement entropy. Entanglement entropy of von Neu- 
mann measures the degree of bipartite entanglement, and 
thus depends on partitions of the system. Multipar- 
tite entanglement entropy thus seems necessary in order 
to characterize total entanglement entropy of a system; 
however, the method is difficult to applied on general 
systems [18]. We suggest one reasonable way to quan- 
tify the total entanglement entropy of a quantum system 
is to sum up von Neumann entanglement entropies of all 
bi-partite subsystems. Since we only concern the entropy 
change in the evolution, we can subtract a ground value 
to regularize any potential infinity. More on regulariza- 
tion will be discussed in later sections. 

To track entropy evolution, a specified setting is 
needed, but yet general enough to characterize the es- 
sential features. The early universe we are investigat- 
ing at roughly can be divided into two stages: inflation- 
ary, and radiation dominated eras. From the requirement 
of unitarity, our model must provide smooth transition 
between those stages; in addition, we will work in tra- 
ditional FRW homogeneous background. An candidate 
with minimal ad hoc assumptions for our mathematical 
model is generalized Chaplygin gas (GCG), which can 
realize the transition smoothly in its equation of state. 
Instead of attaching it as real substance, we only regard 
it as an effective description and convenient way to track 
down entropy evolution. The details of our model, in- 



cluding parameter space, are discussed in section (III A I 



The result from this model, however, we believe to be 



quite general; see section ( VI I 



In the construction of entanglement entropy evolution, 
there are two contributions arose from two distinctive 
ways of dividing the universe. Regarding the entire uni- 
verse as unitary quantum mechanical system consist of 
observable and unobservable regions, von Neumann en- 
tropy then measures the quantum entanglement between 
those two parts of universe due to the evolution of back- 
ground metric, and is thus a homogeneous contribution. 
In this unitary setting von Neumann entropy character- 
izes the probability of different histories, since the to- 
tal amount of information is kept constant by Louvillc's 
theorem. The entanglement entropy is then obtained by 
the standard Wheeler-DeWitt formalism, where configu- 
rations of fields represent different histories in the path 
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integral sense. 

The second contribution comes from the entanglement 
between two opposite momentum sectors of the uni- 
verse, and is entitled inhomogeneous entanglement en- 
tropy. In our model of homogeneous background, it at- 
tributes to entanglement of pairwise matter momentum 
modes, which can be computed by standard cosmologi- 
cal perturbation. The total inhomogeneous entanglement 
entropy is then obtained by summing over the entire mo- 
mentum spectrum. 



III. GENERALIZED CHAPLYGIN GAS 
A. The Model 

When applied to cosmology, the Chaplygin gas [19] 
models the change of cosmic content by regulating the 
equation of state of the background fluid instead of 
the form of the potential. It was first suggested by 
Kamenshchik[20 in an attempt to smoothly interpolate 
the de Sitter phase and the radiation dominant era with- 
out ad hoc assumptions. Generalized Chaplygin gas 
models (GCG) were subsequently suggested [2"T] [2"2~] [2"3~] . 
some of which managed to unify dark energy and dark 
matter [24 . Instead of endowing it as a physical sub- 
stance, we will only invoke it as an effective description. 

The fluid density of GCG is given by [15]: 



P : 



A 



B 



x 4(l+7) 



(1) 



where A and B are positive constants. Note that in 
this expression we have implicitly taken a reference scale 
factor at some point during inflation to be so = 1. 
One can also manifest this reference scale by rescaling 



A i-» Aa 







B i — ^ Ba 



4(1+7) 




Consider the parameters /3 = — 1 and 7 = — 2. Ob- 
serve that when taking the limit a — > 0, one recovers 
from Eq.([T]) the fluid density of de Sitter space charac- 
terizing the inflation era; when taking a — > 00, the fluid 
density becomes that of a gas in the radiation-dominated 
era. Thus the two eras can be smoothly patched under 
this effective description. The constraints on the param- 
eters and 7 to realize this scenario will be discussed in 
Section HHJ 

From the conservation of energy one can deduce the 
pressure for the Chaplygin gas as 



P : 



1 + ,8-4(1 + 7) 
3(1 + 7) 



B 



a 4(l+7)' 



(2) 



After eliminating the scale factor a, one can obtain the 
equation of state for generalized Chaplygin gas. 

One can also model GCG with an underlying minimal 
coupling scalar field <f> by matching the energy momen- 
tum with the fluid described above. Denoting the scalar 



potential by V(^), we have: 

P=\<P 2 + V^), p =lj?-V(<t>)- (3) 

The classical expression for the scalar potential as well 
as the equation of motion can be immediately obtained 
by relating a and 0, as in |19j . For this purpose, we only 
need the limiting case for the early-time and the late-time 
approximations, respectively, with the partition roughly 
coincides with the crossing between the inflation and the 
radiation-dominated epochs. These are 



• Early time [a/a < (A/B) 

, , , 1 fl+J, a 
(a) ~ 0o + -\ +- — In — , 
k y 1 + 7 a 

1 + /3 



(4) 



v(4>) 



1 - 



6(1 + 7) 



A 1 /(i+7) a (V«) g-AW-W (5) 



where A = k^J (1 + /3)/(l + 7) > is the coupling 
constant, tfio = <fi(ao), and n 2 = 8ttG. 



Late time [a/a » (A/B) 



V*i. 



v(4>) 



f - In— , 
k ao 

l B l/(l+7) a 4 e -2K(0-^)^ 



(6) 
(7) 



B. Parameter Space of GCG 



In this subsection we will briefly summarize the con- 
straints on the generalized Chaplygin gas ([!]) 

First we require that the inflation takes place before 
radiation dominant era, which amounts to the condition: 



B 



.4 



(8) 



a 4(i+7) ^ 

as a — > 0. Since both A and B are positive, we have 

5= l + /3_4(l+ 7 ) > 0. (9) 

Two more conditions should be imposed on f3 and 7 so 
as to insure that our model does interpolate between an 
early inflationary phase of the type of quintessence and 
a subsequent radiation dominant phase. These are: (i) 
the energy density must induce a period of inflation and 
(ii) the inflation should not cause a super-inflationary 
expansion; i.e. < H. That is, a super-accelerating 
phase of the universe, where the energy density grows as 
the universe expands, should not occur. 

By combining these two ansatz with the above inequal- 
ity ([9]), we arrive at the following additional constraints: 



1 + /3 



1 + /3 < 0, 
1 + 7 < 0, 

2(1 + 7) > 0. 



(10) 
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In particular, the last constraint results from the require- 
ment for an inflationary period a > as a — > 0, which 
directly leads to the constraint on the coupling constant 
of the GCG scalar held 



X < V2k. 



(11) 



It is possible to have inflation at intermediate value of a 
(see Appendix [A]) . However, in this paper we will only 
consider the case of inflation at the very beginning. 



C. Quantizing Generalized Chaplygin Gas 

Before entering into explicit calculations, we like to 
briefly comment on the problem one often encounters 
when implementing the quantization on the generalized 
Chaplygin gas model, such as the one we introduced in 
Section pII A| When promoting a scalar field to a quan- 
tum field, one expects that the scalar potential would re- 
ceive loop corrections. By inspecting on Eq. one can 
quickly recognize that this model is nonrenormalizable by 
simply noticing the scalar coupling A has mass dimension 
— 1. However, since the scalar field tj> serves as the infla- 
ton in the inflationary period, which decoheres entirely 
sometime after the onset of inflation, we do not expect 
the already decohered field to regain quantum fluctua- 
tions in the remaining period of inflation. We can there- 
fore safely invoke the semiclassical approximation to the 
system beyond the point of decoherence, where one only 
has to consider the finite saddle-point contribution, and 
no regularization is required. On the other hand, before 
the onset of decoherence of the inflaton, one does have 
to treat the scalar field <f> that represents the Chaplygin 
gas fully quantum mechanically, thus some sort of regu- 
larization is still needed. 

In our computation of the homogeneous part of the 
entropy, we invoke the minisuperspace approximation, 
which amounts to truncating large degrees of freedom 
and keeping only the homogeneous modes. Therefore the 
minisuperspace approximation serves as our regulariza- 
tion, based on which a finite result is obtained in Section 



IV D For the correct evaluation of the inhomogeneous 
part of the entropy at early times, one does have to im- 
pose a cutoff scale, below which our model is effective. 
We will discuss more on the issue of quantum corrections 
to the inhomogeneous entropy at early times at Appendix 



IV. ENTANGLEMENT ENTROPY: 
HOMOGENEOUS PART 

As explained in section [TTJ the contribution to the 
cosmological entanglement entropy is divisible into two 
parts; in this section we will address the homogeneous 
contribution using the Wheeler-DeWitt (WDW) formal- 
ism and the minisuperspace treatment. The wavefunc- 
tion of the universe depends on two modes: one from 



gravitation, the other from the scalar field of the Chap- 
lygin gas with the spatial dependence suppressed, thus 
the name "homogeneous" . 

We will treat the gravitational degree of freedom using 
the flat FRW metric: 



ds 2 = -N 2 dt 2 + a(t) 2 dx 2 



(12) 



Throughout this paper we will follow the convention k 2 = 
8nG = 6. 

The action of our model is S = S-eh + "SgcGj where 
the Einstein-Hilbert action can be written as 



S E h = - [d 3 xdtN~a- x Pl 



(13) 



where P a = aa/N is the canonical momentum conjugated 
with a. The spatial integral gives a volume Vo, which 
should not concern us and will be dropped hereafter. The 
matter action is 

Sgcg = / d 3 xdtNa 3 (^0 2 - V(cf>)j , (14) 

where V{<f>) is the scalar potential of the Chaplygin gas. 

Varying with the lapse function N, one can obtain the 
Hamiltonian equation H — 0, which can be quantized 
straightforwardly using the canonical prescription. Here- 
after we will stick to the gauge N = 1. 



Semiclassical Calculation of Entanglement 
Entropy 



As mentioned in Section fill C| since the scalar field 4> 
decoheres sometime during inflation, we can substitute 
by a using the equation of motion. Therefore, the matter 
action is replaced essentially by the fluid density, and the 
Hamiltonian equation boils down to 



H = - l -a- x P 2 a +a 3 p(«) = 0, 



(15) 



which, after the canonical quantization P a <-¥ —i-^, be- 



H^= Giav + 



Matter 



* = 0, 



(16) 



where = \l/(a) is the wavefunction of the universe, " 
denotes differential operators, and ©Matter equals essen- 
tially to the generalized Chaplygin gas density Q. The 
equation above is the Wheeler-DeWitt equation with one 
dynamical degree of freedom a. For later convenience, we 
will factor out a a" 1 factor in both terms, making it re- 
semble the Schrodinger equation for a unit mass particle. 

Since we would like to compute the von Neumann en- 
tropy from the matter wavefunction, which can be ex- 
tracted from the Born approximation [25] 



JSi 



(17) 
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where EH is the gravitational part, which, to the lead- 
ing order, satisfies the Hamilton-Jacobi equation 



( dS EH 
\ da 



- 2a 4 p(a) = 0. 



The matter wavefunction tp satisfies the Schrodinger 
equation corresponding to the potential of the scalar 
field: 



V(<Ka)) = V(a) 



(18) 



To simplify the computation, we will invoke the early- 
time and the late-time approximations, i. e., Eq.([5|and 
Eq.([7]), respectively. In particular, the potential has the 
asymptotic forms: 



V(a) 
V(a) 



V (a/aor^ 1 ^ 
Vq (a/a y 4 , 



Early-time 
Latc-timc 



(19) 
(20) 



with Vq and V ' as well as the criterion of division given 
in Section III A We observe that the potential is posi- 
tive. Therefore the stationary ground state wavefunction 
decays exponentially. 

Since this is a semiclassical treatment, we are only al- 
lowed to obtain the wave function through the WKB ap- 
proximation: 



%j;{a) ~ Ce 



-/». JW{x) 



dx 



(21) 



where a; is a point in time before the semiclassical ap- 
proximation breaks down. (cf. Section IV B[ ) 

From the above equation we can compute the von 
Neumann entropy with the density matrix p(ai,a2) = 
\ip{a 1 )\ 2 5{ai - a 2 ), 



AS(a) 



p In pda . 



(22) 



A bar is used to distinguish it with the fluid density. 

The numerical results are shown in Figure [T] for early 
times and Figure [2] for late times, where the late time is 
entirely in the radiation-dominated epoch. 

From the figures one can observe that the evolution 
of the "homogeneous entropy" beyond the gravitational 
decoherence point is roughly increasing followed by de- 
creasing in radiation epoch. This is expected since par- 
ticle horizon only present in inflationary period, thus en- 
tanglement entropy is expected to decrease after some 
time in the transition. 



B. Breakdown of Semiclassical Treatment 

As has been pointed out, the semiclassical treatment 
of Chaplygin gas, as used in most literature, breaks down 
at the very beginning of the universe. Physically, we ex- 
pect that at the onset of inflation the inflaton should 
execute large quantum mechanical fluctuations that seed 
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FIG. 1: The evolution of the homogeneous entanglement en- 
tropy under the early-time approximation, Eq.Q & Eq.|5]|, 
contributed from the entanglement between the observable 
universe and that outside the horizon, obtained from the 



WKB wavefunction of Wheeler-DeWitt equation Eq.(21 1 
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FIG. 2: The evolution of the homogeneous entanglement en- 
tropy under the late-time approximation, Eq. (JgJ & Eq.|7]l, 
contributed from the entanglement between the observable 
universe and that outside the horizon, within the radiation 
dominant era 



the CMB spectrum. One should therefore treat the scalar 
field quantum mechanically before the decoherence of in- 
flation occurs. 

Consider the scalar potential after substituting the 
saddle point relation, cf> — 4>(a), which at early times 
is of the form V(a) - a - { - 1+ ^^ 1+ ^. The validity of the 
semiclassical (WKB) approximation requires the quan- 
tum potential to vary sufficiently slowly. To be more 
precise, it requires the WKB characteristic length 



WKB 



h 

Pel 



to change slowly over the distance, that is, 

("WKB 



da 



« 1, 



(23) 



(24) 



6 



where without loss of generality we have taken the mass 
to be 1. 



When inserting the early time potential in Eq.( 18 ), we 
find 



dlwKB 




d(V-V2) 


da 




da 



1 + 3-2(1+7) 

a 2<-i-,) _ (25) 



However, with the constraints on parameters in Eq.(lO), 
it is easy to see that the right hand side goes to infin- 
ity when a — > 0, thus violating the validity condition 
of WKB semiclassical approximation. A full quantum 
treatment of both a and <p is therefore necessary for early- 
time universe, which we will carry out in the following 
Section HVCl 



We will implement a fully quantum mechanical treat- 
ment at early times via the full Wheeler-DeWitt equa- 
tion: 



O 



Grav 



o 



Matter 



*(o, 



0, 



(26) 



where we have promoted <j> to be an independent op- 
erator. Written explicitly, we have (in the convention 



1 d d 
o— a- 



1 d 2 

2" da" da 2£hj? 
or, with a = In a, 



+ a b V(0) ¥(a,$=0, (27) 



In order to recover the Chaplygin gas density at the 
saddle point, we adopt the scalar potential V{4>) = V^e^ 
from Eq.(6). Thus the WDW equation is just the scalar 
field with an exponential potential in the FRW metric. 
This model is known to be exactly solvable using the 
following change of variables [26] : 



V = 



/2Vo e Za ~^l 2 
~3~ I - (A/6) 2 

/2Vo e 3Q - A0/2 
"3" 1- (A/6) 2 



cosh 9(a, 



sinh#(a, i 



(29) 



sinh 6{a, cf>) + — cosh 6(a, </>) 



(30) 



where 6(a, (j>) = 3<t>-\a/2. The WDW equation Eq.g 
becomes 

d 2 d 2 



1 



(31) 



_ dx 2 dy 2 
which can be easily solved to give 

(32) 



Take a Gaussian wave packet data with the peak at kg 
and width a, we obtain a wave packet solution 



1*1 



(k 2 - 1)3/V 



a 4 y 2 



x exp 



2_2 



\/k' 2 -lx - fay] 



(*g-l)3 + y3<j4 



(33) 



C. Wavefunction at the Very Beginning of Universe sical solution 



A typical profile of | * | 2 is shown in Figure [3] We can see 
that the distribution is peaked at the saddle point with 
a nonzero quantum-mechanical dispersion similar to the 
free field wavepacket in [2"T] . 

The central value kg can be determined from the clas- 



r a (we take <f> = at a = 0) by 



identifying the location of the "ridge" where the distri- 
bution maximizes at any slice of a given a. The result is 
ko = — 6/\/36 — A 2 . Note that in our convention k = 
A < y/2~K ~ 3.46, which is below 6 and therefore is real 
and finite. 



FIG. 3: Distribution of the wavepacket solution (33 1 

metric inflation a 




D. 



Revised Calculation of Early-time 
Entanglement Entropy 



With the wavefunction distribution, Eq.(33l, at hand, 
we can compute the von Neumann entropy S = 
— /dad^l^l 2 In |*| 2 . Note that if our wavepacket has 
perfect peak with respect to <fi then we will recover the 
equation of motion Q by saddle point approximation 
upon the <p integration, and the entropy will be the same 
as in section HV Al 

The Resulting entropy profile is shown in figure [4] (for 
A = 3.3, o = 0.8). 

The general tendency of the entropy evolution is accel- 
erating, which patches well with the result in Figure [T] 
This is expected since we use a well-defined wavepacket. 

In Figure |4j there's an early period where the numer- 
ical computation becomes less stable; however, in the 
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FIG. 4: Early time homogeneous entanglement entropy evo- 
lution with quantum correction 

Entropy S 




range of tolerance one can observe an entropy minimum 
occurs in the early period. 

We give two explanations here: the entropy minimum 
indicates a decreasing of entropy at the onset of inflation, 
where the quantum fields are out of equilibrium, thus 
thermodynamics are not expected to hold exactly. An- 
other possibility is minisuperspace approximation fails to 
hold at early time, which might be tested by adding back 
degrees of freedom. In our next part of inhomogeneous 
entropy we will pursue this phenomena further. 



V. ENTANGLEMENT ENTROPY FROM 
COSMOLOGICAL PERTURBATIONS 

In this section we will address the part of entanglement 
entropy induced by the entanglement between different 
momentum modes, which we call the "inhomogeneous" 
entropy. This contribution can be calculated using the 
standard cosmological perturbation techniques [28] [29] . 
The calculation in this section largely follows [17] . 



Computation of Entropy from Cosmological 
Perturbation 



Let us denote by v the usual Mukhanov-Sasaki vari- 
able. The comoving curvature perturbation is then 
C(x,i) = v(x,t)/z(x,t), with z(x,t) — a^/e/AG and 



e = -din .fiT/dln a 



(34) 



is one of the slow-roll parameters. The conjugate momen- 
tum of Q is denoted as 7r = d^C- The Gaussian random 
state in a single scalar field inflation is then characterized 
by the covariant matrix between two momentum modes 
k and — k, which is related to the density matrix p by 



C = Tr(pVFt) = ( P C P£« 

r Q-K r 7T 1 



(35) 



where the trace is taken over the functional space, and 

^=(Ck,7T_ k ), 

where ir = 2a 3 e£ is the momentum conjugate with £. 



Each component in Eq. ( 35 ) can be identified (to tree 
level) as: 



AttG 



p c (k,t)Hc(k,t)| 2 , 

Re(C(k,i)ftC(M)*), 



P ff (M) = (— mc(MI 



(36) 
(37) 

(38) 



We can then evaluate the entanglement entropy from cos- 
mological perturbation by the equations: 

S = 2 E k [Ok + 1) ln K + 1) - n k In n k )] , (39) 



{n k + kf = dctC k . 



(40) 



The field ^(k, t) satisfies in Fourier space the equation 



(Ft) 2 " 



v = 0, 



(41) 



which governs the evolution of ((k,t). Here F 
{1/{2V3k))X 2 A 1 / 2 ^ 1 +^ is a constant. 



The basic steps to calculate the entropy from Eq.(39) 



are: (i) evolving the background metric by solving Ein- 
stein equation; (ii) solving the comoving curvature per- 
turbation at time t for different modes k according to 
Eq.(41); (iii) summing the entropy contributions from 



all modes which amounts to integrating Eq. ( 39 1 over the 



momentum space k. In the calculation we employ the 
full tree-level scalar potential without early or late time 
approximation, given explicitly in [19j . 
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FIG. 5: Inhomogeneous entanglement entropy from tree-level 
cosmological perturbation 

The numerical results are shown in Figure [5] Observe 
that the tree-level contribution increases as power law 
in the radiation-dominant epoch started with S'inhom ~ 



10 , which compensates almost entirely the exponen- 



tial decrease of the homogeneous part (section IV A) 
of the entropy, which drops below Shom ~ 10 12 a ftcr 
(a/ao) ~ 10 13 . On the other hand, the homogeneous en- 
tropy during early times continues to increase, resulting 
in the monotonic rise of the net entropy in the semiclas- 
sical regime after the decoherence of the inflaton. 

In the early stage of inflation there exists an entropy 
minimum, which confirms the result in Fig [4] of Sec- 
tion |IVD| that even after adding back some of the trun- 
cated degrees of freedom, there must still remain a per- 
sistent entropy minimum. We can therefore expect that 
the entropy minimum may actually persist to higher or- 
der quantum correction, rendering this turning point a 
real physical feature. We include some discussions about 
higher order correction in Appendix [B| 



VI. CONCLUSION 

In this paper we reconsidered inflation as a possible 
solution to the entropy problem on the basis neglected in 
the past literature that the probability of inflation, when 
taking its accompanied horizon into consideration, is not 
necessary small. Indeed our computation suggests that 
the initial entanglement entropy should have a sizable 
value in order to compensate for its subsequent reduction 
at the beginning of inflation (cf. Fig. [5]), which indicates 
that the onset of inflation could be highly probable. 

We explored this notion by investigating the evolution 
of entanglement entropy from the onset of inflation to 
the radiation era using a toy model, which we believe to 
be generic enough to capture the essential feature of en- 
tanglement entropy evolution because most calculations 
we made do not involve the model's artificial transition 
period between the two epochs, while at the same time 
it has the advantage of preserving unitarity instead of a 
discontinuous phase transition, which is disallowed in the 
quantum mechanical framework for the entire universe. 

Our results of quantum entanglement entropy evolu- 
tion are shown in Figure [l] for homogeneous early evo- 
lution after decoherence of inflaton, Figure [2] for homo- 
geneous late-time evolution, Figure [4] for homogeneous 
early evolution near the onset of inflation, and Figure [5] 
for inhomogeneous evolution at the tree-level. 

The net result for the period after decoherence of infla- 
ton (in our case provided by the generalized Chaplygin 
gas at early time) is that the entanglement entropy in- 
creases monotonically, which results in a large amount of 
entropy generation and agrees with our experience with 
thermodynamic laws. 

On the other hand, the situation before decoherence 
is more subtle. In particular, entropy minimum arisen 
from quantum corrections at early times appears in both 
the homogeneous part (Fig. |4| as well as the tree-level 
inhomogeneous part (Fig. [5]) previously observed in [17] . 
This minimum may result from nonequilibrium effects 
during the early stage of accelerating expansion. The 



persistence of the entropy minimum after putting back 
some of the degrees of freedom suggests that it is a real 
physical feature and we expect this entropy minimum not 
be completely erased by higher order corrections. 

Our results can be summarized as follows: after the on- 
set of inflation the entropy will soon evolve with the ex- 
panding universe dynamically to a minimum and rapidly 
increase hereafter. We therefore suggest that the special 
initial condition required for the thermodynamic arrow 
can be naturally provided by the inflation itself without 
the need to introduce any ad hoc assumptions. 
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Appendix A: Range of Parameters for Inflation 

In this appendix we will show it is possible in our model 
to realize inflation starting from an expanding but non- 
inflationary stage, thereby the parameter space for in- 
flation can be wider than considered in section IIIIBl al- 



though in this paper we restrict to the case of inflation 
at the start. 

For the inflation to start, the kinetic energy must be 
small in comparison with the potential energy [7] , which 
in our model amounts to the situation where the time 
derivative of the field 



1 



_12 = 1 

T 2 da 



1 



<r ~ IP ~ ZC a -(i+«/a+7) 



(Al) 



as well as its spatial derivative, \k 2 a 2 (f) 2 , be small in 
comparison with the potential 

V~(l- I±£) Ca-^+WiW. (A2) 

The former requires 



1 + > 3(1+7), 



(A3) 



while the later will be the case after a 3> exp(l/r) with 
r = 2 — (1 +/3)/(l +7). Comparing to the requirement 
of early inflation at a — > in Eq.(10), we see that there 



is a finite nonzero range of parameters where the kinetic 
term is suppressed as the universe grows but not neces- 
sarily inflates, thereby paving the way for the system to 
become "potential dominated" and ready for the onset 
of inflation. 
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Appendix B: Higher Order Quantum Corrections to 
Inhomogeneous Entropy 

In this appendix we briefly discuss higher order correc- 
tions to the inhomogeneous entropy. Quantum correc- 
tions will only be important at early times. According to 
Eq.([5]), the scalar field potential has the form 



(Bl) 



where we have absorbed the constant <f> into Vq. Thus in 
comparison to the lowest order calculation, where only 
modes k and —k are entangled: 



S = 



, d 3 k 1 
d W £ 2 



a 3 (CfcC-fc - a- 2 k 2 Ck(-k) , (B2) 



there is an additional entanglement due to the nonlinear 
interaction of Eq.(Bl). This contribution can be par- 



tially accounted for by the noise kernel N k (t,t') and the 
diffusion kernel D k (t,t') [29]: 

1 ddet(C) f a 3 e\ 2 
~ \4arGj 



2 dt 
+P<(t) 



— oo 

t 



dt' [D k G(t,t') - N k G iet {t,t')] 
dt' [N k d t G Iet (t,t') - D k d t G(t,t')} 



(B3) 



with N k and D k obtained through the standard diagram- 
matic rules, where higher-order loop effect enters. 



However, upon expansion of the potential Eq.(Bl I one 



finds infinitely many types of vertexes, where each type 
will have its own coupling due to counter terms: 



V, 



quantum 



(</>) = c + citf) + c 2 (f> 2 + 



(B4) 



This immediately leads to problems when performing 
perturbative loop calculations. Therefore numerical lat- 
tice computation seems desirable to obtain higher order 
corrections to early time inhomogeneous entropy. 
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